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Abstract 

We investigate the consistency relation relating the squeezed limit of the bispectrum 
to the scalar spectral index in single field models of inflation. We give a simple integral 
formula for the bispectrum in the squeezed limit in terms of the free mode mode functions 
of the primordial curvature perturbation, in any Lorentz invariant single field model of 
inflation and without resorting to any approximation, generalizing a recent result obtained 
by Ganc and Komatsu in the case of canonical kinetic terms. We use our result to verify 
the consistency relation in an exactly solvable class of models with a non-trivial speed 
of sound. We then verify the consistency relation at the first non-trivial order in the 
slow-varying approximation in general single field inflation (a known result) and at second 
order in this approximation in canonical single field inflation. 



renaux@ ape . univ-paris7.fr 



2 



Contents 

1. Introduction 

2. A reminder of general single field inflation Q 

3. The squeezed limit of the bispectrum in general single field inflation 6 

3.1. The action for short- wavelength modes in the long- wavelength mode backgrounc7 

3.2. Quantizing ( n and applying the Keldysh-Schwinger formalism |9 



4. Verifying the consistency relation llO 

4.1. An exactly solvable class of models : power-law inflation with a constant speed 



of sound 11 
4.2. Slowly varying general single field inflation |l2 

4.2.1. A warm-up: canonical inflation at leading order in the slow- varying 
approximation |l4 

4.2.2. The case of an arbitrary speed of sound at leading order in the 
slow-varying approximation [L5 

4.2.3. Canonical inflation at next to leading order in the slow- varying 
approximation |l8 



21 



5. Conclusions |20 
Acknowledgments |21 

A. Details on the calculations in subsection 13.21 

B. Power law Dirac-Born-Infeld inflation 

C. Useful properties of Hankel functions |23 

D. Useful integrals |24 

E. Details on the correction terms 



22 



24 



References 



3 



1. Introduction 

Observations of the anisotropies of the Cosmic Microwave Background radiation reveal 
that the primordial fluctuations which seeded them are to a good approximation adiabatic 

- in which case they are solely characterized by a primordial curvature perturbation ( 

- scale-invariant and Gaussian [lj]. However, a small amount of non-Gaussianity is still 
allowed by the data and the information contained in this non-Gaussian component will 
contribute to huge advance in our understanding of the early universe by discriminating 
between otherwise competing models (see for instance 0, 0] for recent reviews and 

for more observational aspects). One of the most important realizations in this respect is 
the identification of a consistency relation between the primordial two-point correlation 
function — the power spectrum — and a particular geometrical limit of the three-point 
correlation function — the bispectrum — that is valid in any single field model of inflation 

lim (CbCfaCta) = -(2,r)W 3 )(^k i )(n,(A 1 ) - l)P ( (h)P c (k 3 ) (1.1) 



where 
and 



(CkiCk 2 ) = (27r) 3 5( 3 )(k 1 + k 2 )P ( (h) (1.2) 

is the scalar spectral index. Originally derived by Maldacena in his study of the bispectrum 
generated by a phase of slow-roll single field inflation the relation (11. ip has later been 
generalized by Creminelli and Zaldarriaga to any single field model upon using very general 
arguments [8|. To understand its theoretical and observational relevance, let us remind 
that the observational constraints on the primordial bispectrum are most often quoted in 
terms of the dimensionless momentum- dependent function /tvl(^i, k 2 , k 3 ) defined by 

(C kl Ck 2 Cka) = M^^k^/jvL ( P ^ P ^ + 2 permutations) (1.4) 

i 

Comparing (II. ip and (II. 4p . and using the fact that P((k) oc k~ 3+ ( ns ~ l \ one finds that 

f N 9 L (h) = ^(l-n s (k 1 )) (1.5) 

where we have defined 

f S NL( k i) = Van fNL(h, k 2 , k 3 ) (1.6) 

and where sq qualifies the "squeezed" configuration of momenta under consideration 
(k 3 — > 0,ki ~ k 2 )- Given that the deviation of the primordial power spectrum from 
scale invariance is tightly constrained, n s = 0.963 ± 0.012 (68%CL) [1], any convincing 
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detection of a large bispectrum signal in the squeezed limit (f^ L > 1) would hence rule 
out all models of inflation based on a single scalar field. 

Because of this far-reaching implication, it is very important to understand every aspect 
of the above consistency relation and to verify it by explicit calculations of the primordial 
bispectrum. This has been done for instance in references at first order in a slow- 

varying approximation. Recently however, Ganc and Komatsu gave an integral formula 
for the squeezed bispectrum in single field inflationary models with canonical kinetic terms 



that is valid without any approximation 12J . Although they were not able to derive the 



consistency relation from it, they used their approach to verify the relation (11 .11) non- 
perturbatively in an exactly solvable class of models. 

In this paper, we extend their work in several directions. After reviewing certain aspects 
of the most general Lorentz- invariant single field models of inflation in section [2j we gener- 
alize the main result of Ganc and Komatsu by giving an integral formula for the squeezed 
bispectrum in this general class of models without resorting to any approximation (section 
[3]). We then use this result to verify Maldacena's consistency relation in specific cases in 
section HI We first consider an exactly solvable class of models with a non-trivial speed of 
sound. Second, we verify the consistency relation in general single field inflation at first 
order in the slow-varying approximation. Although this result has already been obtained 



in references [10, lU|, we believe our different derivation to be useful. We demonstrate 
this by finally deriving the consistency relation in canonical single field inflation at second 
order in the slow- varying approximation. We give our conclusions in section [5] and leave 
the details of some long calculations to the appendices. 



2. A reminder of general single field inflation 



In this section, we wet up our notation and we give the second- and third-order scalar 
action in general single-field inflation that were calculated in [10, |l3|] and that we will use 
in the following. 

The action we consider takes the form (we use units in which h = c = M v \ = 1) 



S 



d 4 xy/—g 



-R+P(X, 



(2.1) 



where is the inflaton, X = —\g^ u d^(j)d u (l) and P(X,(p) is the most general Lorentz- 
invariant Lagrangian that is a function of <fi and of its first derivative (see 14Hl7l| for 
extensions to multifield inflationary models). In a spatially flat Friedmann-Lemaitre- 
Robert son- Walker spacetime, with metric 



ds 2 = -dt 2 + a 2 (t)d^ 2 , 



(2.2) 
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where t is cosmic time, the scalar field is homogenous and its energy-momentum tensor 
reduces to that of a perfect fluid with energy density 

p = 2XP X - P (2.3) 

and pressure P. The dynamics of the scale factor and of the inflation field are then 
governed by the Friedmann equations 

3H 2 = p (2.4) 
p=-3H(p + P) (2.5) 

(the equation of motion of the scalar field reducing to the continuity equation (12.51) ). 

In the following, we assume that the scalar field Lagrangian P(X, <fi) is such that a 
prolonged stage of inflation occurs and we study the cosmological perturbations generated 
in such scenarios. In this respect, it is useful to introduce the parameters 

e ~ H*~ 2H* ' V - He' S ~ Hc s ' 1 ] 



E = XP X + 2X 2 P jXX = ^ , {21 



X = X 2 P tXX + ^X 3 P, XXX (2.9) 

where c 2 s is known as the "speed of sound" squared of perturbations, that we require to 
be comprise between and 1 to avoid any pathological behaviour. When the scalar field 
Lagrangian is canonical, P = X — V(4>), the speed of sound equals one and the parameters 
s and A identically vanish. 

To compute the action at second and cubic order in the perturbations, it is useful to 



work in the ADM formalism 18 in which the metric is written in the form 



ds 2 = -N 2 dt 2 + h ij (dx i + N i dt)(dx j + N j dt) . (2.10) 

The lapse iV and the shift N l appear indeed as Lagrange multipliers in the action (12. ip 
and hence can be algebraically expressed in terms of the true physical degrees of freedom. 
Restricting to scalar perturbations, there is only one such quantity, namely, the gauge 
invariant scalar perturbation £ that is constant outside the horizon and of which we want 
to determine the statistical properties. ( is most easily defined in the comoving gauge 
where the inflaton is homogeneous and where the three-dimensional metric hij takes the 
form 

hij = a 2 e^bij . (2.11) 
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Solving the constraint equations for N and N l in terms of ( and plugging them back into 
the action ( 12. ip . one then finds the second-order action 13 



So 



dtd 3 x 



a^C 2 -ae{dCf 



and third-order action 10 



(2.12) 



where 



dtd 3 x 



-a 3 (S(l - i) + 2A)^ + °^(e - 3 + 3c 2 )CC 2 



+-^(e-2s + l-ci)adCY-a 



2c? V c? 



9 2 denotes the inverse Laplacian and, in the last term 

8L . 5S 2 d / , e A 

and 



(2.13) 



(2.14) 



(2.15) 



/(C) 



4c? 



ciH 



Aa?H 2 1 



-(flOCflO + ^WCftCSiC))] 



2a 2 F 



(2.16) 



3. The squeezed limit of the bispectrum in general single field 
inflation 



In this section, we give an explicit integral form for the bispectrum generated during a 



12 



phase of general single field inflation in the squeezed limit, generalizing the result of 
We closely follow their derivation, to which we refer the reader for more details. Let us 
nonetheless now outline the strategy and give important precisions before mowing to the 
calculation itself. 



Our final goal is to compute the bispectrum (Cki (£)Ck 2 (^)Ck 3 (£) ) in the extreme squeezed 
limit in which k 3 — > (and hence k\ ~ k 2 ) at a time t where all the three have become 
classical and have reached their constant value, typically a few-efolds after the largest 
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momenta have crossed the sound horizon at ki 2 c s ~ aH. We will first calculate (CkiCk 2 )k 3 
(the expectation value of Ck x Ck 2 given that Ck 3 has a particular value) and then correlate 
this result with £ ka to find (CkiCk 2 Ck 3 ) = ((CkiCk 2 )k 3 Ck 3 )- For that purpose, we split ( into 
a large-scale, classical, background part (i and a small-scale quantum part ( s : 

0-/ T^U^, ts= [ T^Cke^ (3.1) 



k<k 



(where fc* is chosen such that k 3 < <C k\ ~ k 2 ), so that £ = (i+( s - Introducing the latter 
equation into the second- and third-order action (12.121) and (12 . 131) respectively, we keep the 
terms of order (coming from the second-order action) and (^Q (coming from the third- 
order action). The terms of order £| provide the equations of motion for ( s that determine 
the free mode functions while the terms of order (^Q are treated as perturbations that 
enable one to compute (CkiCk 2 )k 3 using the Keldysh-Schwinger formalism (other terms in 
the third-order action like ( s (f do not contribute to this correlation function). 

We now want to stress an important point : the consistency relation (II. ip relates the 
bispectrum to the scalar spectral index only for very squeezed triangles, i.e. when fc 3 — > 
(see the discussions in 0, 0] for instance). In this limit, the quantum to classical transition 
for the long wavelength mode is pushed away to the past infinity. Therefore, even in cases 
where large interactions occur while the short-wavelength modes are under the horizon 



19l . |20j, such as in models with features, the long- wavelength mode can be treated as 
classical and constant. Hence, although the Keldysh-Schwinger formalism involves an 
integral from early time to the time i at which we evaluate the bispectrum, one can 
neglect, in the limit k% — > 0, the interactions terms containing time derivatives of the long 
wavelength mode, as well as terms containing its spatial gradients obviously 2 . This will 
enable us to drop a number of terms. 

3.1. The action for short- wavelength modes in the long- wavelength mode 
background 

Inserting C = + Cs i n to the second-order action (I2.12p . the second-order action for the 
short-wavelength part ( s is straightforwardly derived: 

S S 2 = I dtd 3 x [a 3 4c/ - ae(d( s ) 2 } . (3.2) 



2 One should be cautious with terms in ( s ( s ,id or (d 2 (^ s )(d 2 ( S: i)(d 2 Q,i), that do arise in our 
calculation, where the spatial and time derivatives act on Q in opposite ways. One can check that these 



terms are indeed negligible in the squeezed limit (the case of the first one is treated in [12j ) . That would 



• 2 

not be the case for instance for terms in £ s d~ 2 Q, absent in our calculation. 



Similarly, at zeroth order in the squeezed limit, i.e. neglecting any time and space deriva- 
tive of (i, the terms of order ( 2 Q m ( [SB]) are found to be 



S s 2 i 



dtd x 



are 



e - 3 + 3c 2 s )dC s 



(If 



a e ( 7] 



+-( e -2s + i- cj)ad( s y + -j- h> occ 
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The last terms in g is most efficiently treated by using a field redefinition 



C — Cn + -5-GCn + • • • 



(3.3) 



(3.4) 



where we have omitted the term with a time derivative of the short wavelength mode since 
this field redefinition will only be evaluated at the time i where all the modes have become 
constant. The second-order action for the redefined field ( n , that we call So, and the cubic 
interaction action between Q and ( n , that we call Si nti (3), then take the form 



S 

Sint,(3) 



6 * 2 

dtd 3 x [a 3 — ( n - ae(d( n ) 2 } 



dtd 3 x 
3 
2 



a 3 e 

IF 



6 - 3 + 3c2)C,C„ + ^(e - 2s + 1 - c 2 Md( n 



are I T] 



(3.5) 



Note that the vertex in ( 3 in (I2.13p . absent in models with standard kinetic terms, and 
which generates equilateral type non-Gaussianities, is manifestly irrelevant in the squeezed 
limit. The non-canonical structure of the action actually only manifests itself in the 
coefficients of the vertex in (I3.5P being different from the canonical case. 

Note that we are eventually interested in the two-point correlation function of ( s , and 
not of ( n . From (13.41) . the link between the two is found to be (this is worked out in details 



in 



12 for Co 



(C, kl (t)C,k 2 (*)> - (Cn, kl (t)Cn,k 2 (t)) +2|^P C (/c 1 )0, kl+k2 



(3.6) 



In the next subsection, we will use the cubic action (13.51) to compute (Cn,ki(*)Cn,k a (^)) i n 
the extreme squeezed limit. 
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3.2. Quantizing £ n and applying the Keldysh-Schwinger formalism 

We now follow the standard procedure to proceed to the quantification of ( n (see 
22]). We first expand ( n in Fourier space 



21 



or 



C " (X) = / (0 C - k6 " X (3 - 7) 
and promote £n,k to a quantum operator 

(n,k(t) = u k (t)a k + ul(t)a[_ k (3.8) 

where and are annihilation and creation operators that satisfy the canonical com- 
mutation rules 

[a k ,4,] = (2n) 3 6^(k-k'). (3.9) 

The so-called mode functions u k (t) (as well as u* k (t)) satisfy the classical equations of 
motion derived from (I3.5P 

o?—^Uk J + aek 2 u k = . (3.10) 



Of 



Equivalently, introducing the canonically normalized field in conformal time r = J dt/a(t) 

v k = zu k , z = , (3.11) 

c s 

equation ( I3.10P can be recast in the familiar form (a prime denotes a derivative with 
respect to r) 

v' k + (c 2 s k 2 - ^jv k = (3.12) 
where v k has to satisfy the quantization (Wronskian) condition 

v l v 'k ~ v kV k = -i (3.13) 

and the appropriate vacuum condition (we will treat the case of the standard Bunch-Davies 
vacuum in the following section). 

From (I3.8P and (13.91) . one deduces the power spectrum of ( n (which is the same as the 
power spectrum of £): 

(Ck(t)Ck'(t)} = (27r) 3 5( 3 )(k + k')K(t)| 2 (3.14) 

so that 

P ( (k)= lim \u k \ 2 . (3.15) 

aH^S>kc s 
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We now turn to the two point correlation function of ( n induced by the cubic interactions 
in (E3D 



(Cn,ki(f)Cn,k 2 (^)) 



0. 



(3.16) 



when kx and — k 2 are different (although very close). At tree- level in the Keldysh- 
Schwinger formalism 23J, |24| , this is given by [25] 



(Cn,ki(*)Cn,k 2 (*))Ci 



k 3 



—I 



dt{0\Cn^(t)Cn,^H m (t)\0} 



C.C. 



(3.17) 



where H 



1,(3) 



-L 



int,(3) 



is the cubic order interaction Hamiltonian and all fields are in the 



interaction picture (which means that they are free fields). Inserting the expression (13. 5p 
into (13.171) . one finds (see appendix El for the details of the calculation): 



with 



™l (r) 



dr 



(Cn,ki(^)Cn,k 2 (^))Ci,k3 - ^0,ki+k 2 



^ e ' - ,2\„2/„./* \2 , 2e 



(3.18) 



=_( e _ 3 + 3^)a a «) a + ^(1 - ci + e - 2 S )a 2 ^«) 2 



+ 



2e (V 



,3„,/* * 



+ C.C. . 



(3.19) 



Using this result together with f)3.6p . and correlating with as announced, one finds 



hm(C kl Ck 2 Ck 3 ) = (2n) 3 5%^) P dh)(Pdh)^ + F) 



cm 



(3.20) 



Equations ( 13 . 19[) and (I3.20j) provide an expression for the bispectrum in the (extreme) 
squeezed limit in terms of the mode functions of ( in general Lorentz-invariant single field 
inflation, the case of a canonical scalar field studied in 12J being recovered when c s = 1 
and s = 0. We stress that the only approximation that we used in deriving them is to 
work at zeroth order in the squeezed limit. The various parameters e, r\ . . . entering into 
(I3.19p . (I3.20p are only short-hand notations and are therefore neither necessarily small nor 
slowly varying. The vacuum state is also left arbitrary. 



4. Verifying the consistency relation 

In the previous section, we derived an expression (13.191) . (I3.20p for the squeezed limit 
of the bispectrum in general single field inflation. We now use this result to explicitly 
verify the consistency relation (II .ip in various cases. We begin by considering an exactly 
solvable class of models in subsection 14.11 We then move to the main calculations of this 
paper: the verification of Maldacena's consistency relation at first order in the slow- varying 
approximation in general single field inflation and at second order in the slow-varying 
approximation in canonical single field inflation. 
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4.1. An exactly solvable class of models : power-law inflation with a con- 
stant speed of sound 

In this subsection, we consider power-law inflation with a constant, but otherwise 
arbitrary, speed of sound, i.e. an inflationary phase with parameters e (< 1 in order 
to realize inflation) and c s constant (and hence r] = s = 0). In appendix [Bj we give 
an example of a model, considered in 0, 27], that realizes such a scenario with a 
Dirac-Born-Infeld type of Lagrangian. Note however that the proof of the consistency 
relation below does not depend on the details of such a realization. 

From e = constant, one deduces that a oc (t — to) 1 ^ an d hence 

( T ) = (-r)"^ (4.1) 
up to a choice of normalization and origin of time. Equation 13.121 then takes the form 

v'l + (cy - ^) v k = (4.2) 



r 2 



with 



v = j — . (4.3) 
2 1-e V ; 



The exact solutions of equation ( 14.2[) . ( 14. 3h are known, and those with positive frequency 
modes — i.e. we choose the Bunch-Davies vacuum — that obey the normalization condi- 
tion (I3.13P read 

Mr) = i^-V^HPi-kcr) (4.4) 

or equivalently 

Mr) = ^ C ^^f e (-TyHW(-kc s T) (4.5) 

where H„ is the Hankel function of the first kind and of order v (useful properties of the 
Hankel functions are collected in the appendix [C]) . From the late time behaviour of the 
Hankel function ( 1C.5j) . one finds 

lim n M^) = -dH>) (4 - 6) 

2v2e7r \kc s J 



and hence the scalar spectral index 

n _ x _ dln(fc 3 |n fc (r)| 2 ) 
dink 



= 3-2u = — . (4.7) 
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Now using the relation (10.3 j) . one calculates that 



Plugging (14.11) . (14. 5 p and (14. 8 p into our general result (13.191) . one then finds 



(4- 



in 



dxx \Hfl x {x)\ + (l-c 2 s + e) / dxx(H {2 \x)) 

-kic a f J —k\c a f 

(4.9) 

Using formulas given in appendix O, one gets 

f°° 2 4 
i I dxx(H^\x)) + c.c. = v , 

Jo 71 



+C.C. 



(4.10) 



so that 



4(21/- 1) + f i - 1 J (3- 21/) . 



p dh) C 2V" - ■ Vc 2 v- (4-11) 

With i/ given in (14.31) . one sees that the terms with negative powers of c s disappear and 
one is left with 

wrr^i- (412) 

It is then clear from equations (13.201) and (14.71) that we have verified the consistency 
relation (II. ip in this exactly solvable class of models with an arbitrary speed of sound, 
thereby generalizing the results obtained in 1_2| for canonical power law inflation. 



4.2. Slowly varying general single field inflation 

In the remainder of this paper, we consider a slowly-varying inflationary phase (though 
with an arbitrary speed of sound unless otherwise specified) with perturbations in the 
standard Bunch-Davies vacuum. The consistency relation (II. ip has already been checked 
in this general class of models at leading order in the slow- varying approximation in 



references 10] and 11] . Although our calculation obviously shares common features with 
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the ones in these papers, we stress that it is largely different. Indeed, Chen et al. 
calculated the full bispectrum through a complete quantum calculation, and later took 
the squeezed limit, whereas we consider this limit directly, resulting in a considerable 
simplification. This latter approach is also followed by Cheung et al. but in the framework 



of the effective field theory of inflation 28(, and in practice with a completely different 



method. Moreover, we will show that our approach enables one to verify quite readily 
Maldacena's consistency relation at second order in the slow-varying approximation in 
models with canonical kinetic terms, a new result to the best of our knowledge. 

By a slowly- varying inflationary phase, we mean that the parameters e, 7] and s, that 
were arbitrary for the moment, are considered both as much smaller than unity, which 
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we note as 0(e), and slowly varying, i.e. jj (e,?), s) = 0(e 2 ). In the following, we refer 
to an expression as being of the n-th order (in the slow-varying approximation) when it 
is accurate up to 0(e n ) terms. As for the speed of sound, we remind the reader that 
although it equals one in canonical single field inflation, there are known examples, like 



DBI inflation [29|, |30| , where it can be much less than unity. In the following, we hence 
leave it arbitrary. 

For our purpose, we will need the solution to the mode equation 13.121 (and that verifies 
the normalization condition (13. 1 3[) ) up to first order in the slow- varying approximation. 
This is given by 3 (we refer the reader to [l(| for an explicit derivation) 

u k {y) = ^_^^(i + £ + |) ,3/^(1) ((1 + e + s)y) (1 + o{e*)) (4.13) 
where we chose the Bunch-Davies vacuum and we defined 

and 

,-§ + « + 2 + i. (4.15) 

Note that at zeroth order in the slow- varying approximation, in which a = — -J^(l + (9(e)) 
and v = | + 0(e), expression (14 . 1 3 j) simplifies to the well known result 

u k {r) = , = = (1 + ikc sK r)e- lkc ^ T th order (4.16) 
^/Ae K c sK k A 

u' k (r) = . = k 2 c 2 sK re- ikCsKT th order (4.17) 

y/AexCsKk 3 

where all parameters are considered as constant at this order and we have chosen to 
evaluate them at the arbitrary pivot point (though near r k for consistency) for K = 
axk. Here and in the following, a subscript K indicates that the corresponding quantity 
is evaluated at sound horizon crossing Kc s k = clkHk- 

From the full result (I4.13p . one deduces the asymptotic value of the mode function up 



to 0(e) order [10 



3 The expression given here differs from the one in [3] by a phase factor which can be considered as 
constant at this order and which therefore is irrelevant. 
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where 



C = 7 -2 + ln2 



-0.73. 



(4.19) 



and 7 = 0.577 • • • is the Euler constant. With our convention, the mode functions become 
real at late times so that 



«2(o) = l<4(o)| = p c (k) 



(4.20) 



df(t)\ ti 
dink 



H k 1 (1 — e k ) 1 f(t k ), one obtains the scalar spectral index at 



and hence, with 
second order in the slow varying approximation: 

n s (k) - 1 = -2e k - rj k - s k - 2e\ - rj k e k - s k e k 

- 2(C + l) Vk e k - C-^ Vk - (C + 2)-^s k + 0{e 



Finally, inserting (I4.20p into the expression (I3.19p . the latter becomes 
F = F 1 + F 2 + F 3 with 

-i I drgi(r)a 2 (u' k \) 2 



F 1 = AP^k^Re 
F 2 = 4P c (h)Re 
F 3 = 4P c (k 1 )Re 



dTg 2 (T)a 2 kl(ul^ 2 

oo 
-oo 



where 



9i(r) = 
92(t) = 



-(3 - 3c 2 s - e) 



'1-cf + e-2s) 



(4.21) 

(4.22) 
(4.23) 

(4.24) 

(4.25) 

(4.26) 
(4.27) 

(4.28) 



In the following, when evaluating (14.231) . (I4.24p and (14.251) . we simply note k = ki for 
brevity. We note also that F 3 is 0(e) smaller than Fi and F 2 and hence can be neglected 
at leading order in the slow-varying approximation. 

4-2.1. A warm-up: canonical inflation at leading order in the slow- 
varying approximation 

As a warm-up, we consider the case of a canonical inflationary Lagrangian, for which 
c s — 1 and s — 0, at leading order in the slow-varying approximation. Using (14. 161) and 
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(I4.17p . one finds, at zeroth order in the slow- varying approximation 



P c (k) 
P c (k) 



—e K kRe 



-^Re 

k 



—i 



2ikr 



dre 



%(l-ikr) 2 e 2ikT 



(4.29) 



As usual in such calculations, the integrals that are required are regulated in the infinite 
past by using the appropriate contour in the complex plane r — > — (oo — it) and their 
integration bound is extrapolated to r = as most of their contributions comes from the 
period around horizon crossing [25| (we will treat a case where this is not true in subsection 
I4.2.3p . Using the integrals given in appendix [D] (where we collect the integrals that are 
needed here and in the following), one finds 



F 1 + F 2 = 2e K P c {k) 



(4.30) 



Remember also that the time f is chosen such that the mode functions have then reached 
their asymptotic value, typically a few efolds after sound horizon crossing. Hence, at 
leading order in the slow- varying approximation, the parameter r](f) appearing in the 
general result (13. 201) can be considered as being equal to t\k (this will not remain true at 
next to leading order as we will see). From (14.301) and ( 13.201) . we thus find 



lim (C kl Ck 2 Ck 3 ) 



-(2nf5^(J2 k i) P dh)Pdh) {~2e K - rj K + 0(e 2 )) 

i 

-(27r) 3 5^(J2 k *) P dki)Pdh) (-2e kl + Vkl -(D(e 2 )) . (4.31) 



It is then clear from the expression of the scalar spectral index (I4.2ip that the consistency 
relation ( 11.11) is indeed verified. 



4-2.2. The case of an arbitrary speed of sound at leading order in the 
slow-varying approximation 

We now move on to the case where the speed of sound is arbitrary, for which we will see 
that verifying the consistency relation at the first non-trivial order requires much more 
work. Indeed, given the forms of the coefficients gi and g 2 in (I4.26P and (I4.27P and of the 
solutions (14.161) and (14.171) . one expects from the previous analysis the right-hand side of 
the relation (13.201) to have contributions of order 0(-^) and 0(4-) (not including the factor 

2a c s 

P^(ki)P^(ks)), whereas the scalar spectral index ( 14.2 ip starts at order 0{e). To check the 
consistency relation at this order, one therefore has to verify that the terms evolving 
negative powers of c s disappear as well as to take into account all the (9(e) corrections to 
this naive reasoning. 
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Working for the moment at leading order in the slow-varying approximation requires 
the same type of integrals as in the previous subsection, the only difference being the 
coefficients: 



= [3(1-0 



- -^(l -c^ + e-2s)^-Re 



c s k 



o 

/ | dre 

oo 



I 



2ikc 3 r 



dr 

— (1 - ikc s T) e 



2 J2ikc s r 



2e - 3s 



(4.32) 



where, from now on, we omit the subscript K when the context is clear. Hence we have 
verified that the leading order terms in O(^) disappear. However, the calculation is 
not consistent at this stage because by treating all the slow-varying parameters in the 
integrands as constant, we neglected 0(e) corrections, which, multiplied by the l/c 2 s — 1 
factor in gi and g 2) compete with the result (I4.32p . 

To proceed further, we thus need to expand the various functions in (I4.23p . (I4.24p 
around the time of sound horizon crossing for our pivot scale. However, the spectral index 
being at least of order 0(e), the precise value of this pivot scale should be irrelevant if 
one wants to to check the consistency relation up to order 0(e). The fact that ax will 
disappear from our final result therefore provides a useful check of our calculation. 



There are three types of 0(e) corrections to the integrands in (14.231) and (14.241) that 
need to be taken into account. They were already given in JjJ and we list them below 
while referring the reader to [lo| for explicit derivations. 



Corrections to the scale factor: 
a(r) = 



Ht/tk) + 0(e 2 ) 



H K T H K T H K T 

Corrections to the coupling constants: 

g(r) = g(t K ) + ^(t-t K ) + 0(e 2 g) 



dg 1 t 
= 9{t k ) - — — In — 
dt Hk tk 

where g collectively stands for gi and g<i . 



0( t 



(4.33) 



(4.34) 



• Correction to the mode functions: 
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Denoting as Auk(r) and Au' k (r) the 0(e) corrections to respectively (I4.16P and 
(Q7j) . we obtain from 



Aul(r) 



2^c^k^ 



x 



(e + s)(x — i) + isx 2 + ( — ( e + ^ + ~ — ^c 2 sj l n 



r 



and 



K 



kc sK e 



x 



-(e + — H — )(1 — -) — iex + sx 2 + {it + -rj — -is — sx)x In — 
2, 2, x 2 2 "^fc 



where x = —kc S KT and we defined 

/i(x) = 



dHl l \x) 
dv 



(4.35) 



(explicit expressions of h*(x) and its derivative in terms of special functions are given 
in the appendix [E]) . 



We then consider all three types of corrections to the two integrations (I4.23P and (14.2 
We give the details of these long calculations in appendix [E] and simply give the final result 
here: 



77(f) + F 



VK -^naive AF 

c 2 sK Pdk) Pdk) 



;i + 0(e)) 



Vk 



'sK 



+ 



2e K - 3s K 



- 1 (2e K + VK - 3s K ) + As K (1 + 0(c)) 



= 2e K + t] K + sk + C{e 2 ) 

= 2e k + Vk + s k + 0(e 2 ) (4.36) 

where we added the term coming from the field redefinition in (I3.2UI) . F na i ve / P^(k) was 
already given in (14.321) and AF/P^(k) denotes the contribution induced by the corrections 
that we listed above. From the expression of the scalar spectral index (14.211) . it is thus 
clear that we have verified the consistency relation 

lim (C kl Ck 2 Cka) = -(2vr) 3 5( 3 )(^k l )K(A ;i ) - l^fc^fo) (4.37) 

i 

up to (9(e) order. 
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4-2.3. Canonical inflation at next to leading order in the slow-varying 
approximation 

We now show that one can verify Maldacena's consistency relation at second order in the 
slow-varying approximation in standard single field inflation with surprisingly few efforts 
given our previous work. First, one needs to compute the 0(e) corrections to the result 
( I4.30P for F± + F 2 , which can be straightforwardly derived from the calculations of the last 
subsection. Second, one has to evaluate F 3 (I4.25p . but being suppressed by 0(e) terms 
with respect to F\ and F2, at leading order only. Finally, one has to take into account the 
fact that the term coming from the field redefinition in (13.201) is evaluated at the late time 
f, and not at horizon crossing. We now treat these three calculations successively. 

• Corrections to F x + F 2 (14301) : 

For our calculation in subsection 14.2.11 evaluating Fi + F 2 at zeroth-order in the 
slow- varying approximation was sufficient to obtain the required result. At next to 
leading order however, one must take into account all the 0(e) corrections to the 
naive behaviour of the scale factor, the mode functions and the coupling constants. 
This is exactly what we have done in the last subsection in the general single field 
case!, with the only difference that <?i(t) = -%(3 — 3c 2 — e) and g 2 (r) = — |Yl — 
(? s + e — 2s) now take their "canonical" value when c s = 1, i.e. — e 2 . Hence, no new 
calculations are required and the results can simply be deduced from the integrals 
in appendix |E] by changing their multiplicative coefficients: the corrections to F% 
(respectively to F 2 ) coming from the scale factor and from the mode functions are 
obtained from the results (lE.lj) and ( TE.8j) (respectively ( TE.9j) and (IE. 16[) ) by making 
the replacement l/c 2 s — 1— > — e/3 (respectively 1/cj. — 1 — > e). As for the correction 
to Fi (respectively to F 2 ) coming from the time variation of the coupling constant, 
it can be deduced from ( IE. 31) (respectively from ( IE. 11 1) ) by making the replacement 
—3 — lj (rj— 4s)+6s — > 2er] (respectively — lj (77— 2s)— 2s — > 2erj). Summing 
all these contributions, one finds 



coupling constant (I4.28P and the mode functions (I4.16p . (14.171) into (I4.25P for c s — 1, 




(4.38) 



• Calculation of F 3 : 

Inserting the expressions at leading order for the scale factor a 




one finds 




(4.39) 
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Here, one can not simply extrapolate the integration bound to r = because the 
integral that is required does not converge. We instead have to keep the f dependence 
of the integral: 

„ K \ m (_™«m +a{ _ 2kf) ) (440) 



P ( (k) \H K r) K ) \ kr 



where 



Ci(x) = -J dt^P- (x>0) 

= 7 + \nx + 0(x 2 ) . (4.41) 

Hence, one obtains 



F, 



3 



V kC+ (-P^-) Vk H-kf) (4.42) 



P c (k) \H K rj K J \H K 7] K J 

up to terms that are negligible for f — > and where the numerical factor C was 
already defined in f)4.19p . 

• Corrections to the naive result from the field redefinition: 

Using (14 .34 p . one relates 7](f) to rjK at next to leading order in the slow- varying 
approximation: 

77(f) = V k(i- -^H-Kf) + 0(e 2 )\ . (4.43) 
V H K rj K J 

Adding all the above contributions, one finds that the terms in ln(— kf ) disappear - 
as expected as the bispectrum does not depend on the time at which it is evaluated after 
all the modes have become constant — and one finds: 

V{f) + pTT) = 2tK + r]K + 2e K + e ^ K ( -1 + 2 (t ~ ln ^f)) 



+ (j^^j Vk (C-lna K ) + G(e 3 ) . (4.44) 



Now using that 



VK = VK>[1 + -^—\n^ + 0(e 2 ) ) (4.45) 

e K = e K > (l + m^— + 0(e 2 )) , (4.46) 

it is manifest that the right hand side of f !4.44p is independent of the pivot scale K at this 
order, which provides a useful check of our calculation which would not have been possible 
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with the two "natural" choices K = k and K = 2k. In the former case, the slow- varying 
parameters appear directly evaluated at the scale of interest at the cost of having a more 
intricate calculation of the integrals while in the latter case, the time integrals are more 
straightforwardly performed and taylor expansions such as 1 14. 46 p are needed in the end. 
Now expressing the right hand side of 1I4.44P in terms of parameters evaluated at the scale 
k, one finds 

V(f) + -^y = 2e k + Vk + 2e 2 k + e K rj K (3 + 2C) + {^^j Vk C + 0(e 3 ) . (4.47) 

Hence it is clear from this expression, the general result 1 13. 201) and the scalar spectral index 
114. 2ip that we have successfully verified the consistency relation fll.ip at second order in 
the slow-varying approximation. 



5. Conclusions 

In this paper, we have generalized a recent result from Ganc and Komatsu [12| by giving 
an explicit formula for the squeezed limit of the primordial bispectrum in any Lorentz- 
invariant model of single field inflation 113. 19[) . 1 13 .20p . We stress that besides working in the 
extreme squeezed limit k$ — > 0, no approximation was used, in particular of slowly- varying 
type, in the derivation of this formula which requires an integral over time involving the 
free mode functions of £. We then used this formula to verify the consistency relation f ll.ip 
in various specific cases, exemplifying that no approximation was made by considering an 
exactly solvable class of models, namely power-law inflation with a constant speed of 
sound. 

Specifying then to the slow-varying regime, we were able to verify the consistency 
relation at the first non-trivial order in general single field inflation. Although this result 



has already been obtained in references [10|, we believe our derivation to be useful. 

is 



One of the advantage of the approach presented here, shared by the proof in [11 
that the interaction in ( 3 in the cubic action 112.131) . dominant for generic configurations 
of the momenta, is irrelevant in the squeezed limit to every order in the slow-varying 
approximation, and hence can be neglected from the very beginning. We have also checked 
explicitly that the pivot scale necessary to verify the consistency relation at order 0(e) 
can be chosen arbitrarily. 

We have finally verified the consistency relation at second order in the slow-varying 
approximation in canonical single field inflation, again considering an arbitrary pivot scale, 
which we think to provide a useful consistency check in such calculations. One should 
also note a subtlety that arises when going beyond leading order in the slow-varying 
approximation: the verification of the r independence of the final result for the bispectrum 
113. 19p . 1 13. 20p . We have indeed shown that one must take into account the fact that the 
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terms coming from the field redefinition is evaluated at late time, and note at precisely 
sound horizon crossing, to obtain the correct, meaningful, result. 

Let us finally note that in multiple field models, the squeezed limit of the primordial 
bispectrum is not related in general to the deviation of the curvature perturbation from 
scale invariance and hence is not necessarily small. If a large bispectrum is detected in the 
squeezed limit, multiple field models — together with the large scale nonlinear evolution 
that they offer [3 lH34J] — will have to be considered seriously. 
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as well as 



Appendix A: Details on the calculations in subsection 13.21 



Inserting the Fourier decomposition (13. 7p into the cubic action (I3.5P , the expression (I3.17P 
becomes 



(Cn,ki(f)Cn,k2(*)}Ci,: 



a 3 e 



ae 



+ 



r- 

a 3 e / r] 



(2tt) 6 

3 + 3c£) (0 1 C kl (t) Ck 2 (*) C qi (*) Cc 2 (*) I o> 

2s + 1 - c 2 s )( qi • q 2 ) (0|C kl (t)Ck 2 (t)C qi (t)C q2 (t)|0) 



(0|Ck 1 (t)Ck a (*)Cq 1 (*)Cq 2 (*)|0> 



+ c.c. 



(A.l) 



where we omit the subscript n on ( n from now on if the context is clear. Using the de- 
composition into creation and annihilation operators (I3.8P together with the commutation 
rules (13.91) . one finds 

i(0|Ck 1 (t)Ck 2 (t)Cq 1 (*)Cq 2 WiO) = i(2nf (^(kx + qi )^(k 2 + q 2 ) + 8® fa + q 2 )c^(k 2 + qi )) 

x [u kl (f)u k2 (t)u* qi (t)u* q2 (t)] , (A.2) 

and similarly when one of the Cqi(^) is replaced by Cqi(^)) an d where we have dropped 
a term proportional to (^^(kx + k 2 ) because h\ ^ — k 2 in our calculation. Assembling 
everything and working at leading order in the squeezed limit, i. e. considering that k x -k 2 ~ 
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-k\ ~ —k%, one gets 

(Cn,ki(£)Cn,k 2 (*))c J ,ic3 = iw fci(*) / ^0,ki+k 2 (*) x 



3 



— (6 - 3 + 3cj) 2«(t)) 2 + -(e - 2s + 1 - c{)2k{ (u* kl (t)Y 



+ c.c. (A.3) 



In this kind of integrals that are familiar in the Keldysh-Schwinger formalism, most of 
the contribution usually come from the period around sound horizon crossing kc s ~ aH 
because the rapid oscillations of the mode functions prior to that epoch usually average 
out. In this case, it is obvious that one can put the term C/,ki+k 2 out of the integral 
when k 3 <C k\ ~ k 2 because d has then reached its final constant value since a long time. 
What are less straightforward are the situations where the time integral in f ]A.3j) has a 



non-zero contribution from the epoch where the short wavelength modes are under the 
sound horizon. This typically arises when there are features in the background evolution 
so that the factor ^(^\ changes suddenly (in less than a Hubble time). In that case, 
0,ki+k 2 ma y wen be still in its quantum regime around the time of the feature and a full 
quantum calculation is required for finite k^ (as well as it is not legitimate to drop terms in 
(i in the interaction action (13.31) ). However, because we are interested in this calculation 
in the limit where the triangle formed by the three wavevectors is infinitely squeezed, i. e. 
&3 — > 0, one can legitimately consider Q as constant in the integral. One then arrives at 
the expressions <KTM . (13T20|) . 



Appendix B: Power law Dirac-Born-Infeld inflation 



We give here an example, considered in 26|, |27| , of a scalar field Lagrangian as well as a 



background evolution for which e (12. 7\i and c s (12.61) are constant. For that purpose, we 
consider a class of Lagrangian inspired by string theory and known as of Dirac-Born-Infeld 
type: 



P 



(yi - 2f{4>)X - l) - V{(j>) . (B.l) 



When derived from string theory, cf) in fIB.lj) represents the location of a D3-brane extended 



along the four usual spacetime dimensions in a special radial direction of a throat amongst 
the six extra dimensions of string theory while f(<f)) is known as the warp factor of the 
throat. Here however, we simply regard (IB. II) as a phenomenological model. 

The most salient feature of the Lagrangian (IB. II) is that its non canonical structure 
imposes a speed limit on the inflaton cf 2 < l/f(<f)). Working out the derivatives in (12. 6ft . 
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one actually finds that the speed of sound c s is given by 



1 - f{4>)p (B.2) 
so that it is both positive and less than unity. Using the Hamilton- Jacobi equations [io[ 

i = -2cM) dJ ^ (B.3) 
3H 2 (4>) = — !— ( — I— - 1 ] + V(<p) where (B.4) 



S(<t>) \c(4>) 

c,W=(l + 4/ W ^M)" 1/2 (B.5) 

one can easily show that requiring e (< 1 to realize inflation) and c s to be constant imposes 
the following specific form for the warp factor and potential 

= (1 - <£) exp (tJ^A (B.6) 

as well as the following background evolution 



TH V2^exp[± x l^-<f)) (B.8) 



2c 



s 



H = H exp[±\l^-<P) . (B.9) 



Appendix C: Useful properties of Hankel functions 

Here, we collect the useful properties of the Hankel functions that we used in the calcula- 
tions of section HJ 

«*) = C - 1 ) 

H ( ?\x) = H^>{x) (C.2) 



d 

dx 



.i —H^\x) + vHpXx) = xH^\{x) , (C.3) 
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dxx (H^\x)) 



2 X 

~ ~2 



2 r 



1 



+ 1 V2 



T(v) (2\ v 

small a; 



7T \ .T 




2 e i(*-f"-|) 



large |z| , |argz| < 7r 



Appendix D: Useful integrals 



o 

- / / dre iT = -1 

-oo 







drre 1 



dre tT ln(-r) = 7 



dre lT r ln(-r) = 1-7 



Re 



00 

-1 I — (1 - rrje 



—00 



-1 



Re 



rfr 

-/ I — (1 - ir) ln(-r)e' J 



00 



-1+7 



C.4) 
C.5) 
C.6) 



D.l) 
D.2) 
D.3) 
D.4) 
D.5) 
D.6) 



Appendix E: Details on the correction terms 

In this appendix, we give the details of the correction terms that are needed in subsection 
12^21 to calculate F (13J9|) up to order 0(e). 

We begin by considering the contribution to Fi coming from the correction to the 
leading-order scale factor a f» — -g-. Plugging the correction terms of f!4.33j) into ( 14 . 2 3 [) 
and evaluating gi and the mode functions to leading order, one finds 



A a F i: 



P c (k) 



6e 



1 kc.Re 



—i 



dr(l — ln(—aKkc s T))e 



2ikc s 



1 a 



1 aK \ 

7 — in ) . 

' 2 ; 



(E.l) 
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Similarly, plugging the correction to g\ (c.f. (14.34)) ) 



into (I4.23[) gives 
• A 9 F i: 



e / 1 \ . , ea \ , r 

^u- 1 ( ^- 4s)+6 ^ 



(E.2) 



A 3 Fi 
P c (*) 



-3(tj -4a) 



- 1 + 6a 



ZcCoRe 



r 

-i / <irln(-axA;c s r)e 2jfcCsT 
•/ —00 



(l_ 1 ) (7 _ ln? )(4 + 63 ) + 3 3 ( 7 - ln? ) 



(E.3) 



As for the corrections to the derivative of the mode function (14.35)) . this contributes 
• A u Fi. 



(» [ - - 1 1 lie 



dr 



T 



e 2ikc 3 T x 



I) 3s, \ T 

+ I tie H ) — ax xm hie , 

1 v 2 2 ; J t k v 2 2 



(e+ 2 + 2 )( x - 1 )-^ + sa;2 



<ix 



1 

2 Vcf 



- 1 



I 7 — In 



ax 



2,(l + 7-lnf)-2a( 7 -lnf) 



1 (2e + 77 + s)Re 



L^o 



dx 



e^fx-z) 



X dx e - ix d h 



x 



dx 



where we have isolated the contributions coming from the first and the last term in 
brackets in the right-hand side of the first equality. Using the asymptotic behaviour 



h*(x) (2 - 7)i - % ln(2x) + £>(x 2 



(E.4) 



one indeed sees that the divergences of the two integrals in the last line of (1E.8|) 
cancel each other. Using 



= (1 - -x - -)e- ix - 7xe"Ei(-2zx) 
ax 2 x 



and making use of the integral 



,00 1 
i J dx Ei(— 2ix) = — , 



(E.5) 



(E.6) 
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one evaluates 



Re 



+00 



dx 



e- 2ix (x-i) - 



+ °° dx ■ d 

-~ lx h*(x) 



x dx 



1 
2 



(E.7) 



so that 



A M F X 3 / 1_ 
P c (k) 2 1 ' 



r/(-l+7-lnf)+2 e ( 7 -lnf)- S (l + 2( 7 -lnf 

(E.8) 



The same procedure for the F 2 term f)4.24p gives 
• A a F 2 : 



A a F 2 



2r I — - 1 ) — R<> 

fcC.o 



(— 1 + ln(— Or-A;c s t))(1 — ikc s r) e 



2 „2ifcc s r 



e --1 1-3 7 -In 



2 / 



(E.9) 



A 9 F 2 : 



A<? 2 (r) 



e --1 ){r ] -2s)-2es 



ln^ 



(E.10) 



A g F 2 
P ( (k) 



(r ] -2s)(--l)-2s 



( v -2s)[--l 



c s k 



Re 



"° dr 

— / 1 — ln(— % fcc s r)(l — ikc s r) 2 e 



2 „2ikc a r 



—00 



2s ](_ 2+ 3 (7 _ lnf) ) 



(E.11) 



A U F 2 : 



A M F 2 
P ( (k) 



kc 



■Re 



dr 



(1 — ikc s r)e 



2ikc a r 



X 



(e + — i) + isx 2 + ( — (e + — + ~)(x — i) — ix 2 s) In — 

V 2 2 / tk 



+(e+| + |)e^*(x) 



i(l- 1 )[ £ ( 2 - 6 ( 7 - 1 nf)) + ,(4-3( 7 -.nf)) 



+ — - 1 ) (2e + 77 + s)Re 



ar 



e- ix (l + ix)/i*(x) 



(E.12) 
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where it is clear from the asymptotic behaviour ( lE.4p that the last integral is con- 
vergent. Using ( IE. 61) . 



7T 



]>*(.r) 2it -" - -(1 + %x)e~ lx - ie lx {\ - ix)E\{-2ix) 



(E.13) 



and 



dx 
x 2 



one finds 

so that eventually 
A U F 2 1 ( 1 



7T 



Re 



2x)cos2x + (2-7ra;)sin2x-si(2a;) =0, (E.14) 



-^e- ix {l + tx)h*{x) 



(E.15) 



P c (k) 2 U 1 



2e2-37-ln^f +r75-37-ln^f +s 



a K 



(E.16) 



Summing the six contributions above, one finds that the terms in In ^ disappear, as 
expected, and one gets 



AF 



- 1 (2e + 7] - 3s) + 4s . 



(E.17) 
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